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Abstract. We review the construction of induced representations of the group 
G = SLjtIR). Firstly we note that G-action on the homogeneous space G/H, 
where H is any one-dimensional subgroup of SL2 (K], is a linear-fractional trans- 
formation on hypercomplex numbers. Thus we investigate various hypercom- 
plex characters of subgroups H. Finally we give examples of induced representa- 
tions of SLj (K) on spaces of hypercomplex valued functions, which are unitary 
in some sense. 
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1. The Group SLsfM) and Its Subgroups 

Let SLjfM) be the group of 2 x 2 matrices with real entries and of determinant 
one [26]. This is the smallest semisimple Lie group. Any matrix in SLjlK] admits 
a (unique) decomposition of the form [ , Exer. 1.14]: 

I ^ h\ _ f (X ^ \ I ^ I cos 4) sin 4) 



c dj VO a / \0 lyV— sine}) cos4) 
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for some values a e (0,oo), v G (—00,00) and c}) £ (— tt, tt]. The formula (1) 
rewritten in a way SLg (M) = ANK is knovm as Iwasawa decomposition [26, § III.l] 
and can be generalised to any semisimple Lie group. 

Each out of the three types of matrices in the right-hand side of (1) forms a one- 
parameter subgroup A, N and K. They are obtained by the exponentiation of the 
respective zero-trace matrices: 

(2) A = {(^; e-)=exp(j _°,) , t e k} , 

(3) N . |(j ;)=oxp(o ;),teM}, 

... ( f cost sint\ t\ _^ , ,1 

- U-sint costj-^'^PU oj'^^(-^'^]j- 

The following simple result have an instructive proof. 

Proposition 1. Any continuous one-parameter subgroup o/SL2(M) is conjugate to one 
of subgroups A, N or K. 

Proof. Any one-parameter subgroup is obtained through the exponentiation 
(5) e*^ = ^ 

of an element X of the Lie algebra 5(2 of SL2 (M). Such X is a 2 x 2 matrix with the 
zero trace. The behaviour of the Taylor expansion (5) depends from properties of 
powers X^ . This can be classified by a straightforward calculation: 

Lemma 2. The square X^ of a traceless matrix X = is the identity matrix 

times + be = — dct X. The factor can be negative, zero or positive, which corresponds 
to the three different types of the Taylor expansion (5) ofe^^. 

It is a simple exercise in the Gauss elimination to see that through the matrix 
similarity we can obtain from X a generator 

• of the subgroup K if (— det X) < 0; 

• of the subgroup N if (— det X) = 0; 

• of the subgroup A if (— det X) > 0. 

The determinant is invariant under the similarity, thus these cases are distinct. □ 

Example 3. The following two subgroups are conjugated to A and N respectively: 

2. Action of SL2(IR) as a Source of Hypercomplex Numbers 

Let H be a subgroup of a group G. Let 11 = G/H be the corresponding homoge- 
neous space and s : O — > G be a smooth section [ , § 13.2], which is a left inverse 
to the natural projection p : G O. The choice of s is inessential in the sense that 
by a smooth map 11 ^ 11 we can always reduce one to another. 

Any g G G has a unique decomposition of the form g = s(a))h, where o) = 
p(g) S £1 and h € H. Note that £1 is a left homogeneous space with the G-action 
defined in terms of p and s as follows: 

(8) g : o) g • cu =p(g * s(cu)), 
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where * is the multiphcation on G . This is also illustrated by the following com- 
mutative diagram: 

^ 9* ^ 



a— 

For G = SL2(M), as well as for other semisimple groups, it is common to con- 
sider only the case of H being the maximal compact subgroup K. However in 
this paper we admit H to be any one-dimensional subgroup. Due to the previous 
Proposition it is sufficient to take H ~ K, N' or A'. Then O is a two-dimensional 
manifold and for any choice of H we define [ , Ex. 3.7(a)]: 

(9) s:(u,v)^-^(^JJ [u,v)eR\v>0. 

A direct (or computer algebra [21]) calculation show that: 
Proposition 4. The SL2(]R) action (8) associated to the map s (9) fs; 

(10) L a)' (^'"^^ ^ 



,c dy"' ' ' V (cu-Hd)2-cr(cv)2 ' (cu + d)2 - (j(cv)2 
where cr — —1, and Ifor the subgroups K, N' and A' respectively. 

The expression in (10) does not look very appealing, however an introduction 
of hypercomplex numbers makes it more attractive: 

Proposition 5. Let an imaginary unit i is such that = a, then the SL2(M) action (10) 
becomes: 

/a b\ Qw + b 

(11) , : w i-^. where W — U+ IV, 

^ ' \c dj cw + d 

for all three cases parametrised by a as in Prop. 4. 

The imaginary unit i^ = — 1 is the very well-known case of the complex num- 
bers. The two-dimensional algebra of numbers x + ey with the imaginary unit 
= 1 is known as split-complex, hyperbolic or double numbers [3,11,31], they 
are one of simplest cases of hypernumbers. The parabolic hypercomplex numbers 
also called as dual numbers are of the form x + ey with the imaginary unit^ such 
that £2=0 [4,9,33]. In cases when we need to consider several imaginary units 
simultaneously we use i to denote any of i, e or e. 

Remark 6. The parabolic imaginary unit £ is a close relative to the infinitesimal 
number e from non-standard analysis [ , ]. The former has the property that its 
square is exactly zero, meanwhile the square of the latter is almost zero at its own 
scale. This similarity is exploit in [4] to produce non-standard proofs of the main 
calculus theorems. A similar property allows to obtain classical mechanics from 
the representations of the Heisenberg group [22]. 

Notably the action (11) is a group homomorphism of the group SLalM) into 
transformations of the "upper half-plane" on hypercomplex numbers. Although 
dual and double numbers are algebraically trivial, the respective geometries in 
the spirit of Erlangen programme are refreshingly inspiring [ Id] and provide use- 
ful insights even in the elliptic case [17]. In order to treat divisors of zero, we need 
to consider Mobius transformations (11) of conformally completed plane [9,19]. 



^We use different scripts of the epsilon: e denote hyperbolic imaginary unit and £ — parabolic one. 
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Physical applications of hypercomplex numbers are scattered through classical 
mechanics [33], relativity [3,31], cosmology [7,9] and quantum mechanics [11,22]. 

Now we wish to linearise the action (8) through the induced representations [12, 
§ 13.2; 13, § 3.1]. We define a map r : G — > H associated to p and s from the 
identities: 

(12) r(g) = (s(a)))"^g, where o) = p(g] e a. 

Let X be an irreducible representation of H in a vector space V, then it induces a 
representation of G in the sense of Mackey [12, § 13.2]. This representation has the 
realisation p^^ in the space of V-valued functions by the formula [±2, § 13.2.(7)-(9)]: 

(13) [p^igmcv] =x(r{g-^*s(a,)))f(g-i - a-),. 

where geG, cueH, HeH and r:G— >H, s:£l^G are maps defined above; 
* denotes multiplication on G and • denotes the action (8) of G on £1. 

In our consideration H is always one-dimensional, its irreducible representation 
is always supposed to be a complex valued character. However hypercomplex 
number naturally appear in the SL2{R) action (11), why shall we admit only i^ = 
— 1 to deliver a character then? 

3. Hypercomplex Characters — an Algebraic Approach 

As we already mentioned the typical discussion of induced representations of 
SLjfK] is centred around the case H = K and a complex valued character of K. A 
linear transformation defined by a matrix (4) in K is a rotation of by the angle 
t. After identification M-^ = C this action is given by the multiplication e'*, with 

= —1. The rotation preserve the (elliptic) metric given by: 

(14) x2+y2 ^ (x + iy)(x-iy). 

Therefore the orbits of rotations are circles, any line passing the origin (a "spoke") 
is rotated by the angle t, see Fig. 1(E). 

Introduction of hypercomplex numbers produces the most straightforward adap- 
tation of this result. 




Figure 1 . Rotations of algebraic wheels, i.e. the multiplication 
by e'-*: elliptic (E), trivial parabolic (Pq) and hyperbolic (H). All 
blue orbits are defined by the identity — i^y^ ~ r^. Green 
"spokes" (straight lines from the origin to a point on the orbit) 
are "rotated" from the real axis. 



Proposition 7. The following table show correspondences between three types of algebraic 
characters: 
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Elliptic 


Parabolic 


Hyperbolic 


= -1 


£2=0 


£2 = 1 


w = X + iy 


W = X + £1) 


w — x + ey 


w = X — iy 


w = X — Ey 


w = X — Ey 


e'^ — cos t + i sin t 


e'^^ = 1 + £t 


e*^^ = cosht + £ sinht 


w|g — WW — x^+y^ 


|w|2 = WW = x2 


|w|^ = WW = x2 — y2 


argw — tan^^ ^ 


argw = ^ 


argw = tanh^^ ^ 


unit circle \w\1 — I 


"unit" strip x = ±1 


wnzY hyperbola w^ = 1 



Geometrical action of multiplication by e^* zs drawn on Fig. Ifor all three cases. 

Explicitly parabolic rotations associated with e"^^ acts on dual numbers as fol- 
lows: 

(15) e^" : a+ eb a+ £(ax + b). 

This links the parabolic case with the Galilean group [^^] of symmetries of the 
classic mechanics, with the absolute time disconnected from space. 

The obvious algebraic similarity and the connection to classical kinematic is a 
wide spread justification for the following viewpoint on the parabolic case, cf . [8, 
33]: 

• the parabolic trigonometric functions are trivial: 

(16) cospt = ±l, sinpt = t; 

• the parabolic distance is independent from y if x ^ 0: 

(17) x2 = (x+£y)(x-£y); 

• the polar decomposition of a dual number is defined by [33, App. C(30')]: 

V V 

(18) u + ev = u(l + £ — ), thus |u + £v|=u, argfu+£v] = — ; 

u u 

• the parabolic wheel looks rectangular, see Fig. l(Po). 

Those algebraic analogies are quite explicit and widely accepted as an ultimate 
source for parabolic trigonometry [8, 27, 33]. However we will see shortly that 
there exist geometric motivation and connection with parabolic equation of math- 
ematical physics. 

4. A Parabolic Wheel — a Geometrical Viewpoint 

We make another attempt to describe parabolic rotations. If multiplication (a 
linear transformation) is not sophisticated enough for this we can advance to the 
next level of complexity: linear-fractional. 

Imaginary units do not need to be seen as abstract quantities. As follows from 
Lem. 2 the generators of subgroup K, N and A represent imaginary units of com- 
plex, dual and double numbers respectively Their exponentiation to one-parameter 
subgroups K, N' and A' of SLjfM] produce matrix forms of the Euler identities 
from the fifth row of the table in Prop. 7. 

Thus we attempt to define characters of subgroups K, N' and A' in term of 
geometric action of SL2(M) by Mobius transformations. The action (11) is defined 
on the upper half-plane and to relate it to unitary characters we wish to transfer it 
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to the unit disk. In the elliptic case this is done by the Cayley transform, its action 
on the subgroup K is: 

. , 1 / l — i\ /cost — sint\ fl i\ /e'^ 



2 V-i 1 



sin t cos t 



1 







The diagonal matrix in the right hand side define the Moebius transformation 
which reduces to multiplication by e^'*, i.e. the elliptic rotation. 
A hyperbolic cousin of the Cayley transform is: 



(20) 



1 



1 











/cosht sinht^ 
l^sinht coshty 

similarly produces a Moebius transformation which is the multiplication by e^^^, 
which a unitary (Lorentz) transformation of two-dimensional Minkowski space- 
time. 

In the parabolic case we use the similar pattern and define the Cayley transform 
from the matrix: 

■ 1 -e^ 
^-e 1 

The Cayley transform of matrices (3) from the subgroup N is: 



(21) 



1 -£ 
-£ 1 



1 + £t 




1 - £t 







t 

£t 



This is not far from the previous identities (19) and (20), however, the off-diagonal 
(1, 2)-term destroys harmony. Nevertheless we will continue a unitary parabolic 
rotation to be the Mobius transformation with the matrix (21), which will not be a 
multiplication by a scalar anymore. 




Figure 2. Rotation of geometric wheels: elliptic (E), two para- 
bolic (P and P') and hyperbolic (H). Blue orbits are level lines for 
the respective moduli. Green straight lines join points with the 
same value of argument and are drawn with the constant "angu- 
lar step" in each case. 



Example 8. The parabolic rotations with the upper-triangular matrices from the 
subgroup N becomes: 

(22) (^''l' :-£^t+£(t2-l). 

This coincides with the cyclic rotations defined in [33, § 8]. A comparison with the 
Euler formula seemingly confirms that sinp t = t, but suggests a new expression 
for cosp t: 

cosp t = 1 — , sinp t = t. 
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Therefore the parabolic Pythagoras' identity would be: 

(23) sinp^ t + cospt = 1, 

which nicely fits in between the elliptic and hyperbolic versions: 

sin^ t + cos^ t = 1, sinh^ t — cosh^ t = — 1. 

The identity (23) is also less trivial than the version cosp^ t = 1 from [ ] (see 
also (16), (17)). 

Example 9. There is the second option to define parabolic rotations for the lower- 
triangular matrices from the subgroup N'. The important difference now is: the 
reference point cannot be — e since it is a fixed point (as well as any point on the 
vertical axis). Instead we take e^^, which is an ideal element (a point at infin- 
ity [33, App. C]) since e is a divisor of zero. The proper compactifications by ideal 
elements for all three cases were discussed in [19]. 
We get for the subgroup N ': 



£ t V t2 

A comparison with (22) shows that this form is obtained by the change t 
t^^. The same transformation gives new expressions for parabolic trigonometric 
fimctions. The parabolic "unit circle" (or cycle [15,33]) is defined by the equation 

— y = 1 in both cases, see Fig. 2(P) and (P'). However other orbits are different 
and we will give their description in the next Section. 

Fig. 2 illustrates Mobius actions of matrices (19), (20) and (21) on the respective 
"unit disk", which are images of the upper half-planes under respective Cayley 
transforms [15, § 8]. 

5. Rebuilding Algebraic Structures from Geometry 

We want induced representations to be linear, to this end the inducing character 
shall be linear as well. Rotations in elliptic and hyperbolic cases are given by 
products of complex or double numbers respectively and thus are linear. However 
non-trivial parabolic rotations (22) and (24) (Fig. 2(P) and (P')) are not linear. Can 
we find algebraic operations for dual numbers, which will linearise those Mobius 
transformations? 

It is common in mathematics to "revert a theorem into a definition" and we will 
use this systematically to recover a compatible algebraic structure. 

5.1. Modulus and Argument. In the elliptic and hyperbolic cases orbits of rota- 
tions are points with the constant norm (modulus): either + orx^—y^. In the 
parabolic case we employ this point of view as well: 

Definition 10. Orbits of actions (22) and (24) are contour lines for the following 
functions which we call respective moduli (norms): 

(25) for N : |u+ £v| = u-^ - V, for N ' : |u + £v|' = . 

V + 1 

Remark 11. (1) The expression |(u, v)| = — v represents a parabolic distance 
from (0, ^) to (u, v), see [15, Lem. 8.4], which is in line with the "parabolic 
Pythagoras' identity" (23). 
(2) Modulus for N' expresses the parabolic focal length from (0, —1) to (u, v) 
as described in [15, Lem. 8.5]. 
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The only straight lines preserved by both the parabolic rotations N and N ' are 
vertical lines, thus we will treat them as "spokes" for parabolic "wheels". Elliptic 
spokes in mathematical terms are "points on the complex plane with the same 
argument", thus we again use this for the parabolic definition: 

Definition 12. Parabolic arguments are defined as follows: 

(26) for N : arg(u + ev) — u, for N ' : arg'(u + ev] — —. 

Both Definitions 10 and 12 possess natural properties with respect to parabolic 
rotations: 

Proposition 13. Lef Wt he a parabolic rotation ofw by an angle t in (22) or in (24). Then: 

jwt!'" = |w|"' , arg'" Wt = arg'" w + t, 
zvhere primed versions are used for subgroup N '. 

All proofs in this and the following Sections were performed through symbolic 
calculations on a computer. See Appendices A- B for details. 

Remark 14. Note that in the commonly accepted approach [33, App. C(30')] para- 
bolic modulus and argument are given by expressions (18), which are, in a sense, 
opposite to our agreements. 

5.2. Rotation as Multiplication. We revert again theorems into definitions to as- 
sign multiplication. In fact, we consider parabolic rotations as multiplications by 
unimodular numbers thus we define multiplication through an extension of prop- 
erties from Proposition 13: 

Definition 15. The product of vectors wi and W2 is defined by the following two 
conditions: 

(1) arg''' (W1W2) — arg''' Wi + arg''' W2; 

(2) IW1W2I'" =|wi|"' -Iwal'". 

We also need a special form of parabolic conjugation, which coincides with sign 
reversion of the argument. 

Definition 16. Parabolic conjugation is given by 

(27) u+ ev = -u + £v. 

Obviously we have the properties: |w|''' — |wi''' and arg'''w = — arg'''w. 
A combination of Definitions 10, 12 and 15 uniquely determine expressions for 
products. 

Proposition 17. The parabolic product of vectors is defined by formulae: 

(28) /orN: (u,v) * (u',v') = (u + u', (u + u')^ - (u^ - v)(u'2 - v']]; 

(29) /orN': (u, v) . (u',v') ^ (^ /^ - l) . 

\u + u' (u-hu'J^ J 

Although both expressions look unusual they have many familiar properties: 

Proposition 18. Both products (28) and (29) satisfy the following conditions: 

(1) They are commutative and associative; 

(2) The respective rotations (22) and (24) are given by multiplications with a dual 
number with the unit norm. 

(3) The product W1W2 is invariant under respective rotations (22) and (24). 

(4) For any dual number w the following identity holds: 

|ww| = |wj^ . 
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In particular, the property (3) will be crucial below for an inner product (35), 
which makes induced representations unitary. 

6. Invariant Linear Algebra 

Now we wish to define a linear structure on which would be invariant under 
point multiplication from the previous Subsection (and thus under the parabolic 
rotations, cf. Prop. 18(2)). Multiplication by a real scalar is straightforward (at least 
for a positive scalar): it should preserve the argument and scale the norm of a 
vector. Thus we have formulae for a > 0: 

(30) Q-(u,v) = (u,av + u2(l - q)) forN, 

(31) Q-(u,v) = (u,^^-l^ forN'. 

On the other hand the addition of vectors can be done in several different ways. 
We present two possibilities: one is tropical and another — exotic. 

6.1. Tropical form. Let us introduce the lexicographic order on M^: 

either u < u'; 



u, v) -< (u , V ) if and only if 



or u = u , V < V 



One can define functions min and max of a pair of points on respectively. Then 
an addition of two vectors can be defined either as their minimum or maximum. 
A similar definition is used in tropical mathematics, also known as Maslov dequan- 
tisation or R^iiu and Rj,jax algebras, see [28] for a comprehensive survey. It is easy 
to check that such an addition is distributive with respect to scalar multiplica- 
tions (30) — (31) and consequently is invariant under parabolic rotations. Although 
it looks promising to investigate this framework we do not study it further for now. 

6.2. Exotic form. Addition of vectors for both subgroups N and N ' can be defined 
by the common rules, where subtle differences are hidden within corresponding 
Definitions 10 (norms) and 12 (arguments). 

Definition 19. Parabolic addition of vectors is defined by the following formulae: 

(,, argf'wi • jwil''' +arg(''w2 • Iwsl'" 

(32) arg''(wi+W2) = ■ jj] , 

|wi + 

(33) IW1+W2I'" = |wi|"'±|w2l"', 
where primed versions are used for the subgroup N '. 

The rule for the norm of sum (33) may look too trivial at the first glance. We 
should say in its defence that it nicely sits in between the elliptic |w + w'| ^ |w| + 
|w'| and hyperbolic [w + w'| ^ |w| + |w'| triangle inequalities for norms. 

The rule (32) for argument of the sum is not arbitrary as well. From the Sine 
Theorem in the EucUdean geometry we can deduce that: 

. , |w|-sin(\|;-i|j') |w'| • sin(xjj - 1];'] 
sin(4) J = ■ ■ , sm(i|; - 4)] = ■ , 



where ij;''' — argw''' and — arg(w + w'''). Using parabolic expression (16) for 
the sine sinp 6 = we obtain the arguments addition formula (32). 

A proper treatment of zeros in denominator of (32) can be achieved through a 
representation of a dual number w = u + ev as a pair of homogeneous polar coor- 
dinates [a, r] — [|w|''' • arg''' w, |w|'''] (dashed version for the subgroup N '). Then 
the above addition is defined component-wise in the homogeneous coordinates: 

wi +W2 = [qi + a2,ri +72], where wi = [ai,ri]. 
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The multiphcation from Defn. 15 is given in the homogeneous polar coordinates 



Thus homogeneous coordinates Hnearise the addition (32)-(33) and multiplication 
by a scalar (30). A transition to other more transparent coordinates shall be treated 
withing birational geometry framework [24]. 

Both formulae (32)-(33) together uniquely define explicit expressions for addi- 
tion of vectors. However those expressions are rather cumbersome and not really 
much needed. Instead we list properties of these operations: 

Proposition 20. Vector additions for subgroups N and N ' defined by (32)-(33) satisfi/ 
the following conditions: 

(1) They are commutative and associative. 

(2) They are distributive for multiplications (28) and (29); consequently: 

(3) They are parabolic rotationally invariant; 

(4) They are distributive in both ways for the scalar multiplications (30) and (31) 
respectively: 



To complete the construction we need to define the zero vector and the inverse. 
The inverse of w has the same argument as w and the opposite norm. 

Proposition 21. (N) The zero vector is (0, 0) and consequently the inverse of (u, v) 
is (u, 2u^ — v). 

(N ') The zero vector is (oo, — 1) and consequently the inverse of (u, v) is (u, — v — 2). 

Thereafter we can check that scalar multiplications by negative reals are given 
by the same identities (30) and (31) as for positive ones. 

Remark 22. The irrelevance of the standard linear structure for parabolic rotations 
manifests itself in many different ways, e.g. in an apparent "non-conformality" of 
lengths from parabolic foci, that is with the parameter cr = in [15, Prop. 5.12.(iii)]. 
An adjustment of notions to the proper framework restores the clear picture. 

The initial definition of conformality [ , Defn. 5.11] considered the usual limit 
■y' — )- ij along a straight line, i.e. "spoke" in terms of Fig. 1. This is justified 
in the elliptic and hyperbolic cases. However in the parabolic setting the proper 
"spokes" are vertical lines, see Fig. 2(P) and (?')/ so the limit should be taken along 
them [15, Prop. 5.13]. 



We discussed above various implementations of hypercomplex unitary char- 
acters. Now we can return to consideration of induced representations. We can 
notice that only the subgroup K requires a complex valued character due to the 
fact of its compactness. For subgroups N ' and A' we can consider characters of 
all three types — elliptic, parabolic and hyperbolic. Moreover a parabolic character 
can be taken either as algebraic (15) or any of two geometric (22) and (24). There- 
fore we have seven essentially different induced representations, which multiply 
types to eleven (counting flavours of parabolic characters). 

Example 23. Consider the subgroup H = K, then we are limited to complex valued 
characters of K only. All of them are of the form Xk^ 



by: 



wi • W2 = [air2 + a2ri , rir2] , 



where Wi — [at, Vi]. 



a ■ (wi + W2 ) = a • Wi + a • W2 



(a + b) • w = a • w + b • w. 



7. Conclusion: Induced Representations 




where k e Z . 
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Then the corresponding induced representation (13) p^, act on complex valued 
functions in the upper half-plane = SLjfKj/K as follows: 

These are unitary representations from the discrete series [ , § IX. 2]. 

Proposition 24. Let fk(w) = (w — i)'^/or k = 2, 3, . . . , then 

(1) fk ;s an eigenvector for any operator p^fh), where h e K, with the eigenvalue 
Xk(h] [26, § IX.2]. 

(2) The function K(z,w) — p^(s(z))fk(w), where s(z) is defined in (9), is the 
Bergman reproducing kernel in the upper half-plane [13, § 3.2]. 

Similarly we can get the Cauchy kernel for the limiting case k = 1 of the mock 
discrete series [26, Ch. IX]. There are many other important connections of repre- 
sentation (34) with complex analysis and operator theory. For example, Mobius 
transformations of operators lead to Riesz-Dimford functional calculus and asso- 
ciated spectrum [14]. 

Example 25. In the case of the subgroup N there is a wider choice of possible 
characters. 

(1) Traditionally only complex valued characters of the subgroup N are con- 
sidered, they are: 

'I 0^ 



where t e 



The corresponding induced representation acts on the space of complex val- 
ued functions on the upper half-plane M'^, which is subset of dual numbers 
as a homogeneous space SL2(M)/N. The corresponding formula is 

Cf irr ^ /. Tcv \ /aw + b\ _^ /a b 

p;r In f(w = cxp 1 f , where w = u + ev, q ~ \ , 

^ ^ V cu-h dy Vcw + d/ ^ \^ ^ 

The mixture of complex and dual numbers in the same expression is con- 
fusing. 

(2) The parabolic character Xt with the algebraic flavour is provided by mul- 
tiplication (15) with the dual number: 

^ ^ e"^ = 1 + ert, where t e M. 

Then the corresponding induced representation is defined on the space of 
dual numbers valued functions on the upper half-plane of dual numbers 
by the formula: 

p,(g)f(w) = l + £ — f ' 



cu + d / \ cw + d , 

where w, t and g are as above. This expression contains only dual numbers 
with their usual algebraic operations. Thus it is linear with respect to them. 
(3) The geometric character Xt is given by the action (22). Then the corre- 
sponding representation acts again on the space of dual numbers valued 
functions on the upper half-plane of dual numbers as follows: 

„, ^ ( 2tcv \ ^/aw + b\ TCV (tcv]^ 
p9(g)f(w)= l + £— — f — +—— + £-' 



cu+d/ \cw + d/ cu+d (cu + d)2' 

where w, t and g are as above. This representation is linear with respect to 
operations (30), (32) and (33). 
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All characters in the previous Example are unitary, the first two in a conven- 
tional sense and the last one in the sense of Prop. 18. Then the general scheme of 
induced representations [12, § 13.2] implies their unitarity in proper senses. 

Theorem 26. All three representations of SL2(]R) from Example 25 are unitary on the 
space of function on the upper half-plane K'^ of dual numbers with the inner product: 



(35) (fi,f2 



^ , , - , , du dv 

fi(wjf2(wj — - — , wherew — u+ zv, 



and we use 

(1) the conjugation and multiplication of functions' values in algebras of complex and 
dual numbers for representations and respectively; 

(2) conjugation (27) and multiplication (28) of functions' values for the representa- 
tion p?. 

c 

The inner product (35) is positive defined for the representation but is not 
for two others. The respective spaces are parabolic cousins of the Krein spaces [1], 
which are hyperbolic in our sense. 

There are many important questions to be investigated for those induced repre- 
sentations: relations with the three main series of representations (discrete, princi- 
pal, complementary) of SL2(M] [2()], connections with various SI2 modules [10,29], 
applications to analytic functions [13] and partial differential equations [25], as- 
sociated functional calculi and spectra [14], etc. These directions can be viewed 
as parts of the Erlangen programme at large [17,23] and shall be considered else- 
where. 
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Appendix A. Output of Symbolic Calculations 

Here are the results of our symbolic calculations. The source code can be ob- 
tained from this paper [20] source at http : //arXiv . org. It uses Clifford algebra 
facilities [ l^] of the GiNaC library [?]. The source code is written in noweb [30] 
literature programming environment. 
Calculations for subgroup N and straight spokes 
Reference point: ( — 1 ) 
Reference point arg: 
Reference point norm: 1 

Cayley of the matrix X: 



1 + e°x Ix 

1 - e°x 



Rotation by x: ( u + x 2ux + v + x^ ) 
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Rotation of (uq, Vq) by x: ( x — 1 + x^ ) 
Parabolic norm: — v 
Parabolic argument: u 

Real number ti as a dual number: ( — ti ) 

Product: ( u' + u u'^ + 2u'u + + v'u^ - u'^u^ - v'v + u'^v ) 

Product by a scalar: ( u u'^ + av — au'^ ) 



V 



Real part: ( — uv — + u'^ 
Imag part: ( 1 1 + uv — u'^ ) 
Zero vector: ( ) 
Negative vector: ( u 2u^ ^ ) 

Dual number from argument and norm n: ( — n ) 
Dual number from argument Qi and norm n: ( Qi af — n ) 
Dual number from argument Qi and norm n — norm: n 
Dual number from argument Qi and norm n — arg: Qi 

LinCOmboftwOVectOrsa*(l,0)+b*(-l,0): ( b--b-^-a--3b^a+cf-2ba-3ba- 

P is the sum 9l(P) and J(P): true 

The real part of a real dual number is itself: true 

norm is invariant under rotations: true 

Product is invariant under rotations: true 

Product WW is norm squared: true 

Product (u,v) * (uqjVo) is (u,v): true 

Add is commutative: true 

Add is associative: true 

S-mult is commutative: true 

S-mult is associative: true 

S-mult is distributive 1: true 

S-mult is distributive 2: true 

Product is symmetric (commutative): true 

Prod is associative: true 

Product is distributive: true 

Calculations for subgroup N ' and straight spokes 
Reference point: ( oo — 1 + oo^ ) 
Reference point arg: 
Reference point norm: 1 

■ l-e°x 
Ix 



1+ux l+2ux+u^x'^ 
1 -1+x' 



Cayley of the matrix X: 

Rotation by x 
Rotation of (uq, vq) by x: 
Parabolic norm: 
Parabolic argument: ^ 
Real number ti as a dual number: ^ oo 

Product: ^^T^^^p 

Product by a scalar: ( u ~— ^-7^— ^ ) 



Real part: 00 — 

Imag part: ( 1 -^^^r^ ) 
Zero vector: (00 — 1 ) 
Negative vector: ( u — 2 — v 
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Dual number from argument and norm n: ( oo 
Dual number from argument qi and norm n: ( 



Dual number from argument Qi and norm n — norm: n 
Dual number from argument Qi and norm n — arg: Qi 

Lincomboftwovectorsa* (l,0)+b* (-1,0): ( ^^^^b-lf"^ 

P is the sum m{P) and 3(P): true 

The real part of a real dual number is itself: true 

norm is invariant under rotations: true 

Product is invariant under rotations: true 

Product WW is norm squared: true 

Product (u,v) * (uo,vo) is (u,v): true 

Add is commutative: true 

Add is associative: true 

S-mult is commutative: true 

S-mult is associative: true 

S-mult is distributive 1: true 

S-mult is distributive 2: true 

Product is symmetric (commutative): true 

Prod is associative: true 

Product is distributive: true 

Elliptic case of induced representations 
map r(M): 



character: 



Moebiusmap: c a v- + c ou- + b cu+ b cl+ d au _^^bcy_^av^ 
r \^ 2cdu+d''+C'^u'' + c''v'= 2cdu+d''+C'^u'= + c''v'' 

Moebius map is given by the imaginary unit: true 

Parabolic (N') case of induced representations 
1 



d^ 


cd 


d2 + c2 


d2 + c2 


cd 


d^ 


d2 + c2 


d^ + c^ 




bc^ + c+bd- 


it 


d^ 




d^ + c^ 
d 



2cdu+d'' + c^u^ c(cTJ-+d)v 

2cdu+d2+c2it2 + c2v2 2cdu+d2+c2u2 + c2v2 
c(cu+d)v (cu+d)^ 



2cdu+d2+c2u2 + c2v2 2cdu+d^ + c2u2 + c2v2 

cav^ + cau^+bcu+bd+dau bcv— dav 



map r(M): 



f 1 



- b 

maps-i(M): ( ^ ^ 



character: 



1 



cv 
cu+d 



Moebius map: ^ 



1 

-gu bcv— dav 



cu+d (cu+d)- 

Moebius map is given by the imaginary unit: true 



Hyperbolic case of induced representations 
map r(M): 



d" cd 

d2-c2 d2-c2 

cd d^ 



\ d^ — c^ d^— c^ / 
/ 1_ bc^ + c-bd^ 

maps-i(M): f J 
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2cdu+d''+c^u^ c(cu+d)v 



character: 



2cdu+d^+c^u^ — c'-'v^ 2cdu+d^ + c^u 
c(cu+d)v (cu+d) 



2cdu+d2+c2u2-c2v2 2cdu+d2H 



Moebius map: ^ 
Moebius map is given by the imaginary unit: true 



cav^ — cau^— bcu— bd— dau bcv— day 

2cdu+d^+c^u^ — c^v^ 2cdu+d^ + c^u^ — c^v^ 



Appendix B. Program for Symbolic Calculations 

This is a documentation for our symboHc calculations supporting this paper. 
You can obtain the program itself from the source files of this paper [ ] at arXiv . org; 
KTjX compilation of it will produces the file parab-rotation.nw in the current 
directory. This is a noweb [30] code of the program. It uses Clifford algebra facili- 
ties [18] of the GiNaC library [2]. 

This piece of software is licensed under GNU General Public License (Version 
3, 29 June 2007) [6]. 

B.l. Calculation and Tests. This Subsection contains code for calculation of vari- 
ous expression. See [ i ] or GiNaCinfo for usage of Clifford algebra functions. 

B.1.1. Calculation of Expressions. Firstly, we output the expression of the Cayley 
transform for a generic element from subgroups N and N'. 

16a (Show expressions 16a) = (23e)16b> 
ex XC=canonicaUze.clifford{{TC*X*TCr).evalm()); 

formiilaj}ut{"Ca.ylej of the matrix X: " , XC.subs{sign=0).normal()); 

Uses f ormula.out 16c 17h 17h 17h 17h 18a 26d. 

Then we calculate Mobius action of those matrix on a point. 

16b (Show expressions 16a) += (23e) < 16a 16c > 

dualjiumber W{cliffordjnoebius_map{XC, P.to^matrixQ, e).subs{sign=0).normal{)), 

Wl = W.subs{lst{u=ul, v=vl)); 
formulajDut("Rota.tion by x: ", W); 
if (not W.isjequal{P.rot{x))) 

COllt<^"*** dualnumber : :rot() gives wrong answer *** \\("<C 
P.rot{x) <C "\\) " <C endl <C endl; 

Uses f ormula_out 16c 17h 17h 17h 17h 18a 26d. 

Next we specialise the above result to the reference point. 

16c (Show expressions 16a) += (23e) < 16b 16d > 

/ormM/fl_OMf("Rotation of \\((u_0, v_0)\\) by \\(x\\): ", 

W.swiis(lst(u = uO, V = vO)).subs{Inf= pow{y, -\)).normalQ.subs{y = O).normal{)); 

Defines: 

f ormula_out, used in chunks 16, 17, and 21. 

The expression for the parabolic norm. 

16d (Show expressions 16a) += (23e) < 16c 17a > 

/07'mM/fl_0Mt("Parabolic norm: " , P.normQ); 
/ormM/fl_DMt("Parabolic argument: ",P.argQy, 

Uses f ormula.out 16c 17h 17h 17h 17h 18a 26d. 
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Embedding of reals into dual numbers. 

17a (Show expressions 16a) += (23e) < 16d 17b > 

possymbol U("tl","t_l"), fll("al", "a_l"), n{"n"); 

formidaj}ut{"Real number \\(t_l\\) as a dual number: ", dual_number(n)); 

Uses f ormula.out 16c 17h 17h 17h 17h 18a 26d. 

The expression for the product of two points. 

17b (Show expressions 16a) += (23e) < 17a 17c > 

formula j)ut{"Pioduct: ",P*P1); 

Uses f ormula.out 16c 17h 17h 17h 17h 18a 26d. 

The expression of the product of a point and a scalar. 

17c (Show expressions 16a) += (23e) < 17b 17d > 

formulaj)ut{"Pioduct by a scalar: ",a*P); 

Uses f ormula.out 16c 17h 17h 17h 17h 18a 26d. 

Expressions for the real and imaginary parts. 

17d (Show expressions 16a) += (23e) <l 17c 17e > 

formulajDut{"Real part: " ,P.realjpartQ); 
formula j}ut{"lmag part: " , P.imag^partQ); 

Uses f ormula.out 16c 17h 17h 17h 17h 18a 26d. 

The expression for a sum of two points is too cumbersome to be printed. 

17e (Show expressions 16a) += (23e) < 17d 17f > 

/ /formula_out("Add is: ", (P+Pl).normal()); 

Uses f ormula.out 16c 17h 17h 17h 17h 18a 26d. 
Expression for zero is 

17f (Show expressions 16a) += (23e) < 17e 17g > 

formula j)ut{" Zero vector: " , zerojduaLnumberQ); 

Uses f ormula.out 16c 17h 17h 17h 17h 18a 26d. 
Expression for zero is 

17g (Show expressions 16a) += (23e) < 17f 17h > 

formulaj)ut{"}}egative vector: ",P.neg()); 

Uses f ormula.out 16c 17h 17h 17h 17h 18a 26d. 
Expression for zero is 

17h (Show expressions 16a) += (23e) < 17g 18a > 

formulaj}ut{"Bual number from argument \\(0\\) and norm \\(n\\): ", 

dn^omjirgjnod(Q, n)); 
dual-number PP=dn^omjirgjnod{al, n); 

formtdaj)ut{"T)ual number from argument \\(a_l\\) and norm \\(n\\): ", 
PP.normalQ); 

formula j)tlt{"'Dual number from argument \\(a_l\\) and norm \\(n\\) norm: ", 

PP .normQ.normalQ); 

formulaj3Ut{"'D'aal number from argument \\(a_l\\) and norm \\(n\\) arg: ", 

PP .argQ.normalQ.normalQ); 



Defines: 

f ormula_out, used in chunks 16, 17, and 21. 
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Linear combination of points (1,0) and (—1,0) with coefficients a and b, for the 
linearisation presented in [21]. 

18a (Show expressions 16a) += (23e) < 17h 

formulaj)ut{"LirL comb of two vectors \\(a*(l, 0)+b*(-l, 0)\\): ", 
(fl*dual_number(l,0)+fo*duaLnumber(-l,0)).wormfl/()); 

Defines: 

f ormula_out, used in chunks 16, 17, and 21. 

B. 1 .2. Checking Algebraic Identities. In this Subsection we verify basic algebraic prop- 
erties of the defined operations. 

A dual number is the sum of its real and imaginary parts. 
18b (Check identities 18b) = (23e) 18c> 

test.out{"\\(P\\) is the sum \\(\\Re(P)\\) and \\(\\Im(P)\\) : ", 
P-(ex_to<dual_number>(P.reflZ_pflrt())+ex_to<dual_number>(P.nwfl^_pflrt()))); 

Defines: 

test.out, used in chunks 18, 19, and 21f. 

A dual number maid out of a real a has the norm of real part equal to a. 

18c (Check identities 18b) += (23e) < 18b 18d > 

testj)ut{"The real part of a real dual number is itself: ", 

ex_fo<dual_number>(dual_number(fl).reflLpflrt()).non;(()-fl); 

Defines: 

test.out, used in chunks 18, 19, and 21f. 

The norm is invariant under parabolic rotations, i.e. they are in agreement with 
Defn. 10. 

18d (Check identities 18b) += (23e) < 18c 18e > 

test. out{" norm is invariant under rotations: " , P.norm{)-W.norm{)); 

Uses test.out 18b 18c 18f 27. 

The product W1W2 is invariant under rotations. Prop. 3. 

18e (Check identities 18b) += (23e) < 18d 18f > 

tesL0Mf(" Product is invariant under rotations: " , P*Pl.conjiigateQ-W*Wl.conjugate{)); 

Uses test.out 18b 18c 18f 27. 

Product WW is (0, |wp),Prop. 4. 

18f (Check identities 18b) += (23e) < 18e 18g > 

fesLoMf("Product \\ (w\\bar{w}\\) is norm squared: ", 
{P*P.coniugate{)-dnJromjirg.mod(ArgO, pow(P.norm(), 2)))); 

Defines: 

test.out, used in chunks 18, 19, and 21f. 

The reference point is unit under multiplication. 

18g (Check identities 18b) += (23e) < 18f 18h > 

test.out{"Pioduct \\((u, v)*(u_0, v_0)\\) is \\((u, v)\\): ",P*PO-P); 

Uses test.out 18b 18c 18f 27. 
Addition is commutative. Prop. 1. 
18h (Check identities 18b) +EE (23e) < 18g 19a > 

test.out{"kdd is commutative: " , (P+P1)-(P1+P)); 



Uses test.out 18b 18c 18f 27. 
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Addition is associative. Prop. 1. 

19a (Check identities 18b) += (23e) < 18h 19b > 

tesLoMfC'Add is associative: " , {{P+P1)+ P2)-(P+(P1+P2))); 

Uses test.out 18b 18c 18f 27. 

Multiplication by a scalar is commutative. 

19b (Check identities 18b} += (23e) < 19a 19c > 

tesL0Mf("S-mult is commutative: ",P*a-a*P); 

Uses test.out 18b 18c 18f 27. 

Multiplication by a scalar is associative. 

19c (Check identities 18b) += (23e) < 19b 19d > 

tesLoMf("S-mult is associative: " , b*P*a-a*P*b); 

Uses test.out 18b 18c 18f 27. 

Distributive law a(wi + W2) ~ awi + aw2. Prop. 4. 

19d (Check identities 18b) += (23e) < 19c 19e > 

tesLoMt("S-mult is distributive 1: " , a*{P+Pl)-{a*P +a*Pl)); 

Uses test.out 18b 18c 18f 27. 

Distributive law (a + b)w = aw + bw. Prop. 4. 

19e (Check identities 18b) += (23e) < 19d 19f > 

tesLoMf("S-mult is distributive 2: " , P*(a+b)-{P*a + P*b)); 

Uses test.out 18b 18c 18f 27. 

Product is commutative. Prop. 1. 

19f (Check identities 18b) += (23e) < 19e 19g > 

tesLoMf("Product is symmetric (commutative): " , P*P1-P1*P); 

Uses test.out 18b 18c 18f 27. 
Product is associative. Prop. 1. 

19g (Check identities 18b) += (23e) < 19f 19h > 

fesL0Mf("Prod is associative: " , {P*P1)*P2-P*(P1*P2)); 

Uses test.out 18b 18c 18f 27. 

Product and addition are distributive. Prop. 2. 

19h (Check identities 18b) += (23e) < 19g 

tesLoMtC'Product is distributive: " , {P+P1)*P2-{P*P2+P1*P2)); 

Uses test.out 18b 18c 18f 27. 



B.2. Induced Representations. Here we calculate the basic formulae for Section 2. 
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B.2.1. Encoded formulae. This routine encodes the map s : ^ SL2(M) (9). 

20a (Induced representations routines 20a) = (22a) 20b > 

ex s_mflp(const ex & u, const ex & c) { 
return matrix(2, 2, lst(i',u,0,l)); 



ex s_mflp(const ex &P) ^ 
if {P.nopsQ = 2) 

return sjnap{P.op{0), P.op{l)); 
cerr <C "sjnapO error: parameter should have two operands" <g; end/; 
return s_map{P,l); 



This routine encodes the map r : SL2(IR) H (12). The first parameter is an 
element of SL2(M), the second — is a generic element of subgroup H. We look spe- 
cific matrix of the form K which makes the product MK belonging to the image of 
sjitapi), i.e. its (2, 1) element should be zero. 
20b (Induced representations routines 20a) += (22a) < 20a 20c > 

ex r_mflp(const ex <& M, const ex&cK) { 
ex Kl=K.evalm{), K2; 

1st vars = (^s_fl<symboI>(^Cl.op(2)) ? lst(Kl.op{2)) : lst{Kl.op(l))); 
if (zs_fl<symbol>(ra.op(3))) I 
vars = vars.append{Kl.op{3)); 

K2 = Kl.siibs(lsolve{lst{(M*Kl).evalm0.op{2)=0), vars)).subs{Kl.op{3)=l); 
y else 

K2 = Kl.subs{lsolve{lst{(M*Kl).evalm0.op{2)=0), vars)); 
return pow{K2, -l).evalm{); 



This is the inverse s ^ of the above map s. 
20c (Induced representations routines 20a) += (22a) <i 20b 20d > 

ex p_mflp(const ex & M, const ex & J<C) { 
ex X = K.evalm().op{3); 
exMK={M*pow{rjnap{M,K),-1-))-evalm{); 
ex D=MK.op{3).subs{x=l).nornial(); 

return matrix(l, 2, lst{{MK.op{l).subs{x=l).nornial{)^D).normal{), 

(MK.op{0).subs(x=l).normal{)~^D).normal{))); 

y 



This is a matrix form of the above inverse map pjnapQ. 
20d (Induced representations routines 20a) += (22a) < 20c 

ex p_mflp_m(const ex & M, const ex & X) { 
return {M*poiv(rjnap{M,K),-l)).evalni{); 

y 



B.2.2. Caculation of induced representation formulae. Firstly we define a generic ele- 
ment M of SL2(M). 

20e (Induced representations 20e)= (24d)20f> 
ex M=matrix(2,2, \st(a,b,c,d)), H; 

We consider the three cases. 

20f (Induced representations 20e) += (24d) < 20e 21a > 

sfnM^CflSes[]={"Elliptic", "Parabolic (\\ CN~\\prime\\) ) ", "Hyperbolic"}; 
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In the those cases subgroups holds a generic element of a subgroup H, see (4), (6) 
and (7). 

21a (Induced representations 20e) += (24d) < 20f 21b > 

ex SMfo§-roups=lst(matrix(2, 2, lst(x,-i/,i/,x)), 
matrix(2, 2, lst(l,0,i/4)), 
matrix(2, 2, lst{x,y,y,x))); 

Now we run a cycle over the three cases. . . 

21b (Induced representations 20e) += (24d) < 21a 21c > 

for(int i=0; i<3; i 
H=subgroups[i]; 

COM t <C cases [;] <C " case of induced representationsWW" ^ ewd/; 
//formula.out("M*H: ", (M*H).evalm()); 

Uses f ormula_out 16c 17h 17h 17h 17h 18a 26d. 

. . . and output expression of r (12),. . . 

21c (Induced representations 20e) += (24d) < 21b 21d > 

formula.out{" map \\(r(M)\\): " , r.map{M,H)); 

Uses f ormula_out 16c 17h 17h 17h 17h 18a 26d. 

. . . matrix form of the inverse s^^ (12),. . . 

21d (Induced representations 20e) += (24d) < 21c 21e > 

formula j}ut{" map \\(s"-[-l}(M)\\) : " , p_map.m{M,H).subs{a=(l+b*c)^d).normalQ); 

Uses f ormula.out 16c 17h 17h 17h 17h 18a 26d. 

. . . expression for the argument of the character in (13),. . . 

21e (Induced representations 20e) += (24d) < 21d 21f > 

formula j[>ut{" character : ", r_map{M*s.map{P),H)); 

Uses f ormula.out 16c 17h 17h 17h 17h 18a 26d. 

. . . and finally the action (8) of SL2 (K) on the homogeneous space. 

21f (Induced representations 20e)+= (24d)<21e 
/ormM/fl_OMf("Moebius map: ", p_map{M*sjnap{P.tojnatrixQ),H)); 

testjout("Moehius map is given by the imaginary unit: " , pjnap(M*sjnap(P),H) - 

Clifford jnoehius jnap{a*one , b*one, none, d*one,P.tojnatrixQ,e).subs{sign=i-l)); 
cout <Si {latexout? "\\vspace{2mm}\\hrule" : 

" g„^;. 

} 

Uses f ormula.out 16c 17h 17h 17h 17h 18a 26d and test.out 18b 18c 18f 27. 

B.3. Program Outline. Here is the outline how we use the above parts. 

B.3.1. Test program outline. Firstly we load dual_number support. 

21g (*21g)= 22a > 

#include <cycle.h> 
#include <dualnum.h> 
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The rest of the program makes all checks. 

22a C21g)+= <21g 
(Definition of variables 22b) 
(Test routine 23a) 

(Induced representations routines 20a) 
(Main procedure 24a) 



B.3.2. Variables. These realsymbols are used in our calculations. 

22b (Definition of variables 22b) = (22a) 22c > 

const numeric half{\,2); 

Defines: 

numeric, used in chunks 25 and 26. 



Variables vs oftenly appear under square roots of the form Vl + 2v. To facilitate 
the simplifications of the type (Vl + 2v)^ — 1 + 2v we (falsely) define them to be 
positive symbols. 

22c (Definition of variables 22b) += (22a) < 22b 22d > 

possymbol f("v"), ul("v' v2{"v' ' "); 

Other real variables. 

22d (Definition of variables 22b) += (22a) <i 22c 22e > 

realsymbol !(("u"), m1("u"'), u2("u' ' "), 
fl("a"), b("b"), c("c"), d("d"), x("x"), i/("y")/ 

Finally this variable keeps the signature of the metric space. 

22e (Definition of variables 22b) += (22a) < lid lit > 

sign("s", "Wsigma"); 

This an index used for the definition of Clifford units. 

22f (Definition of variables 22b) += (22a) < 22e 22g > 

varidx mM(symbol("mu", "\\mu"), 1), nM(symbol("nu", "\\nu"), 2); 

Three generic points which are used in calculations. 

22g (Definition of variables 22b) += (22a) < 22f 22h > 

dualjiumber P{u, v), Pl{ul, vl), P2(u2, v2); 

Here we define a parabolic Clifford units one, eO. 

22h (Definition of variables 22b) += (22a) < 22g 

exe = Clifford Mnit{mu, diagjiiatrix(lst(sign))), 
eO = e.siibs{mu=0), 
one = dirac^ONEQ, 

e2 = cliffordMnit{nu, diagjnatrix(lst{-l,sign))); 
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B.3.3. Test routine. This routine make the same sequence of checks for both cases 

of subgroups N and N'. 

First we define the reference point {uO,vO). 



void parab_rotjub(hool N, bool S) { 

cout <C "Calculations for subgroup \\(N" 

< (N? "" : ""') <C "\\) and " 

< (S ? "straight" : "geodesies") < " spokes\\\\"<ewd/; 
set^ubgroup(N)} 

set straight ^poke{S); 
exX, 

WO=dn^romjirg.mod{ArgO, 1), 

uO=WO.op(0), 

vO=WO.op{l), 

PO=matrix(l, 2, lst(wO, vO)); 

Defines: 

parab_rot_sub, used in chunk 24c. 
These two matrices define the Cayley transform and its inverse. 

(Test routine 23a) += (22a) < 23a 23c > 

coMi < "Reference point: \\(" < PO < "\\)\\\\" < end/; 

cotii < "Reference point arg: \\(" < ex_to<dualjiumber>(WO).flr§-() <C "\\)\\\\" < ewd/; 
cotii < "Reference point norm: \\(" <C exJo<dualjiumber>(WO).norm() < "\\)\\\\" <C end/; 

For the subgroup N we consider upper-triangular matrices, for N' — lower- triangular. 

(Test routine 23a) += (22a) <i 23b 23d > 

if (subgroupJsT) 

X=matrix(2, 2, lst(one, one*x, 0, one)); 

else 

X=matrix(2, 2, lst(oMe, 0, one*x, one)); 
Two different types of Cayley transforms. 

(Test routine 23a) += (22a) < 23c 23e > 

ex TC,TCI; 
if(S){ 

TC=matrix(2, 2, lst(one, -eO, -eO, one)); 
Ta=matrix(2, 2, lst(one, eO, eO, one)); 
} else{ 

TC=matrix(2, 2, lst(one, -eO*half, -eO*tialf, one)); 
TQ=matrix(2, 2, lst(one, eO*half, eO*half, one)); 



(Test routine 23a) 



(22a) 23b > 



> 



Common part of test routine. 



(Test routine 23a) += 
(Show expressions 16a) 
(Check identities 18b) 

cout ^ (latexout? "\\vspace{2imn}\\hrule 



(22a) <23d 



II 



) < endl; 



y 
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B.3.4. Main procedure. It just calls the test routine, calculates the induced represen- 
tation and draws a few pictures. 
We output formulae in KTgX mode. 

24a (Main procedure 24a) = (22a) 24b > 

int mainQi 
latexout=tiue} 

Defines: 

main, never used. 

Preparation of output stream. 

24b (Main procedure 24a) += (22a) < 24a 24c > 

cout <C boolalpha; 
if (latexout) 

emit <g; latex; 
(Drawing pictures 24f) 



Now we call the test routine for both N and N' subgroups. 
24c (Main procedure 24a) += (22a) < 24b 24d > 

parabj-ot^ub{ti\ie, true); 
parabj'ot^ub(fnlse, true); 

/ / parab_rot_sub(false, false); To work with geodesic spokes 
Uses parab_rot_sub 23a. 

Calculation of induced representations formulae. 

24d (Main procedure 24a) += (22a) < 24c 24e > 

(Induced representations 20e) 

And we finishing by drawing several pictures for Figs. 1 and 2. 
24e (Main procedure 24a) += (22a) < 24d 

> 

B.4. Drawing Orbits. To draw cycles we use cycle library [18]. 
Elliptic orbits (circles). 

24f (Drawing pictures 24f) = (24b) 25a > 

ofstream flsi/mpfoie("parab-rot-data. asy"); 
Asymptote <g; "path [] K="; 
for(int i=0; i<6; 

cycle2D(lst(0,0),e2.sufos(sf§-M=-l),!*i*.04) 

.asyjpath{asymptote, -1.5, 1.5, -2, 2, 0, {i>0)); 
asymptote <g; " ; " <S endl; 
asymptote <g; "path [] Kb="; 
cycle2D(lst(0,0),e2.s!(fos(sf^n=-l),l) 

.asy_path{asymptote, -1.5, 1.5, -2, 2, 0, false); 
asymptote < " ; " < endl; 
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Hyperbolic orbits. 

25a (Drawing pictures 24f) += (24b) < 24f 25b > 

flsympfote <C "pathD A="; 
for(int i=0; i<6; { 
cycle2D(lst(0,0),e2.SMbs(s;§7j=l),-!*;*.04) 

.asy.path(asymptote, -1.75, 1.75, -1.5, 2, 0, (!>0)); 
cycle2D(lst(0,0),e2.SMfo(sz^M=l),z*f*.12) 

.asy.path{asymptote, -1.75, 1.75, -1.5, 2, 0, true); 

> 

asymptote <g; " ; " <S endl; 
asymptote <C "path [] Ab="; 
cycle2D(lst(0,0),e2.SMbs(s!§-n=l),-l) 

.asy_path{asymptote, -1.75, 1.75, -1.5, 2, 0, false); 

asymptote <C " ; " <C endl; 

Hyperbolic orbits for reflected orbits. 

25b (Drawing pictures 24f)+= (24b)<25a25c> 
asymptote <g; "path [] At="; 
for(int i=0; i<6} { 

cycle2D(lst(0,0),e2.sufos(sf§-M=l),i*/*.04) 

.asyjpath{asymptote, -1.75, 1.75, -1.5, 2, 0, (/>0)); 
cycle2D(lst(0,0),e2.siite(sz^M=l),-z*i*.12) 

.asyjpat}i{asymptote, -1.75, 1.75, -1.5, 2, 0, true); 

} 

asymptote ^ " ; " <g; end/; 
asymptote <C "path [] Abt="; 
cycle2D(lst(0,0),e2.si(fcs(sf^w=l),l) 

.asyjpath{asymptote, -1.75, 1.75, -1.5, 2, 0, false); 
asymptote <g; " ; " <C end/; 

Parabolic orbits, subgroup N . 

25c (Drawing pictures 24f) -|-= (24b) < 25b 26a > 

asymptote <C "path [] N="; 
for(int /=0; /<6; f-i^) 

cycle2D(l,lst(0,numeric(l,2)),numeric(z,2)-l,e2.SMbs(s/^M=0)) 
.asyjpath{asymptote, -1.5, 1.5, -2, 2, 0, (!>0)); 
asymptote ^ " ; " <C end/; 
asymptote <S "path [] Nb="; 
cycle2D(l,lst(0,numeric(l,2)),-l,e2.su/;s(s/^n=0)) 

.asyjpath{asymptote, -1.5, 1.5, -2, 2, 0, false); 
asymptote ^ " ; " <C end/; 



Uses numeric 22b. 
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Parabolic orbits, subgroup N'. 

26a (Drawing pictures 24f)+= (24b) < 25c 26b > 

flsympfote <C "path[] Nl="; 
for(int i=0; i<5} 

cycle2D(.5*!*/*;+l,lst(0,numeric(l,2)),-l,e2.siite(si^fj=0)) 
.asy-path{asymptote, -1.5, 1.5, -1.5, 2, 0, {i>0)); 
asymptote <C " ; " ^ endl; 
asymptote <§; "path [] Nlb="; 
cycle2D(l,lst(0,numeric(l,2)),-l,e2.sute(s!^M=0)) 

.asyjpath{asymptote, -1.5, 1.5, -1.5, 2, 0, false); 
asymptote ^ " ; " ^ endZ; 



Uses numeric 22b. 

Parabolic orbits, subgroup N geodesic version. 

26b (Drawing pictures 24f) += (24b) < 26a 26c > 

asymptote <^ "pathD Ng="; 
for(int i=0; i<6; 

cycle2D(0.5,lst(0,numeric(l,2)),numeric(!,2)-.5,e2.SMfos(s!gM=0)) 
.asyjpat}i{asymptote, -1.5, 1.5, -2, 2, 0, (z>0)); 
asymptote <C " ; " <S ewd/; 
asymptote <C "path [] Ngb="; 
cycle2D(0.5,lst(0,numeric(l,2)),-.5,e2.SMfos(sfgn=0)) 

.asyjpath{asymptote, -1.5, 1.5, -2, 2, 0, false); 
asymptote <C " ; " <g; ewd/; 



Uses numeric 22b. 

Parabolic orbits, subgroup N' geodesic version. 

26c (Drawing pictures 24f) += (24b) < 26b 

asymptote <S "path [] Nlg="; 
for(int !=0; i<5; i++) 

cycle2D(.25*/*!*!+.5,lst(0,numeric(l,2)),-.5,e2.SMfos(s/^«=0)) 
.asyjpath{asymptote, -1.5, 1.5, -1.5, 2, 0, (/>0)); 
asymptote <g; " ; " <C ewd/; 
asymptote <t; "path [] Nlgb="; 
cycle2D(.5,Ist(0,numeric(l,2)),-.5,e2.SHbs(si^M=0)) 

.asy_path{asymptote, -1.5, 1.5, -1.5, 2, 0, false); 
asymptote ^ " ; " ^ endZ; 

flsymptote.c/ose(); 



Uses numeric 22b. 



B.4.1. Output routines. We use standardised routines to output results of calcula- 
tions. 

26d (Output routines 26d)= 27 > 

\oiA formula j)ut{string S, const ex & F, bool lineend) { 

cout < S < (latexout ? "W (":"")< F<C (latexoiit ? "\\)" : "" ); 
if (lineend) 

cout <C (latexout ? "WW" : "" ) < end/; 

else 

cowt < " ; "; 



Defines: 

f ormula_out, used in chunks 16, 17, and 21. 
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This routine is used to check identities. 

(Output routines 26d)+= 

void testj)ut{string S, const ex & T) ■[ 

cout < S <C (latexout ? "Wtextbf {" : "*" ) 

<C (!s_fl<dual_number>(T) ? ex_to<dual_number>(T).«ormfl/().zs_zero() : 

T.evalmQ.normalQ.isjzerojnatrixQ) ^ {latexout ? "}\\\\" : "*" ) 
< endl; 

y 



Defines: 

test.out, used in chunks 18, 19, and 21f. 
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